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A detailed description and theoretical analysis of experiments achieving coherent coupling between 
an ion Coulomb crystal and an optical cavity field are presented. The various methods used to 
measure the coherent coupling rate between large ion Coulomb crystals in a linear quadrupole 
radiofrequency ion trap and a single-field mode of a moderately high-finesse cavity are described 
in detail. Theoretical models based on a semiclassical approach are applied in assessment of the 
experimental results of [P. F. Herskind et al, Nature Phys. 5, 494 (2009)] and of complementary 
new measurements. Generally, a very good agreement between theory and experiments is obtained. 
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I. INTRODUCTION 

Cavity Quantum Electrodynamics (CQED) consti- 
tutes a fundamental platform for studying the quan- 
tum dynamics of matter systems interacting with elec- 
tromagnetic fields 0, 0. For a single two- level quan- 
tum system interacting with a single mode of the elec- 
tromagnetic field of a resonator, a particularly interesting 
regime of CQED is reached when the rate, g, at which 
single excitations are coherently exchanged between the 
two-level system and the cavity field mode exceeds both 
the decay rate of the two-level system^ 7, and the rate, 
K, at which the cavity field decays pj]. This so-called 
strong coupling regime was investigated first with sin- 
gle atoms in microwave and optical cavities [H, Hj] and 
recently with quantum dots [f| 0] and superconduct- 
ing Josephson junctions [H, In the optical domain, 
the use of ultrahigh-finesse cavities with a very small 
modevolume allows for reaching the confinement of the 
light field required to achieve strong coupling with sin- 
gle neutral atoms 

0, EH- With charged particles, 
however, the insertion of dielectric mirrors in the trap- 
ping region makes it extremely challenging to obtain suf- 
ficiently small cavity modevolumes, due to the associ- 
ated perturbation of the trapping potentials and charg- 
ing effects [ill, [f| . Although the strong coupling regime 
has not yet been reached with ions, single ions in opti- 
cal cavities have been successfully used for, e.g., probing 
the spatial structure of cavity fields [3], enhanced spec- 
troscopy [ll| , the generation of single photons [ljl fl7| , 
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the investigation of cavity sideband cooling [18| . or the 
demonstration of a single ion laser [l9[ . 

For an ensemble of TV identical two-level systems si- 
multaneously interacting with a single mode of the elec- 
tromagnetic field, the coherent coupling rate is enhanced 
by a factor y/N 0. This leads to another interesting 
regime of CQED, the so-called collective strong coupling 
regime where the collective coupling rate g^ = gx^N 
is larger than both k and 7. This regime, first explored 
with Rydberg atoms in microwave cavities [20I ] , has been 
realized in the optical domain with atomic beams 
atoms in magneto-optical traps pll - l^ . Bose-Einstein 
condensates [HI, l27j - and, recently, with ion Coulomb 
crystals [28|. This cavity-enhanced collective interaction 
with an ensemble has many applications within quantum 
optics and quantum information processing [29^, includ- 
ing the establishment of strong nonlinearities I30L 13111 . 
QND measurements [32]-[34|, the production [35l. l36| 
and storage [37], HH of single-photons, t he g eneration of 
squeezed and entangled states of light [2lTl39l liol ] and 



atoms 



2ol 4lj | , t he observation of cavity optomechani- 
cal effects [22ll^. l4^ - |45| . cavity cooling [MEl, and the 
investigation of quantum phase transitions [47] . 

This paper provides a detailed description and a theo- 
retical analysis of experiments achieving collective strong 
coupling with ions |28j. The various methods used to 
measure the coherent coupling rate between large ion 
Coulomb crystals in a linear quadrupole radiofrequency 
ion trap and a single field mode of a moderately high- 
finesse cavity {F ~ 3000) are described in detail. Theo- 
retical models based on a semiclassical approach are ap- 
plied in assessment of the experimental results of Ref. [281 ] 
as well as of complementary new measurements. Gener- 
ally, a very good agreement between the theoretical pre- 
dictions and the experimental results is obtained. As 
also emphasized in Ref. [28|, the realization of collec- 
tive strong coupling with ion crystals is important for 
ion-based CQED (48[ and enables, e.g., for the realiza- 
tion of quantum information processing devices, such as 
high-efficiency, long-lived quantum memories [37l.l49j and 
repeaters [H(| . In addition to the well-established attrac- 
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tive properties of cold, trapped ions for quantum infor- 
mation processing |5ll . l52j . ion Coulomb crystals ben- 
efit from unique properties which can be exploited for 
CQED purposes. First, their uniform density under lin- 
ear quadrupole trapping conditions [53l - l55j makes it pos- 
sible to couple the same ensemble equally to different 
transverse cavity modes (56j and opens for the realization 
of multimode quantum light- matter interfaces foTj , where 
the spatial degrees of freedom of light can be exploited 
in addition to the traditional polarization and frequency 
encodings (58l - [60j . Second, their cold, solid- like nature 
combined with their strong optical response to radiation 
pressure forces and their tunable mode spectrum (6f"l - [6l| 
make ion Coulomb crystals a unique medium to inves- 
tigate cavity optomechanical effects [64]. Ion Coulomb 
crystals could, for instance, be used as a model sys- 
tem to study the back action of the cavity light field on 
the collective motion of mesoscopic objects at the quan- 
tum limit, as was recently demonstrated with ultracold 
atoms 0, HH, |47j . In addition, novel classical and quan- 
tum phase transitions could be investigated using cold 
ion Coulomb crystals in optical cavities [65l - [68j . 

The paper is organized as follows: Sec. [TT| presents 
the theoretical basis for the CQED interaction of ion 
Coulomb crystals and an optical cavity field. The cav- 
ity field reflectivity spectra, and the effective number of 
ions interacting with the cavity field are derived and the 
effect of temperature on the collective coupling rate is dis- 
cussed. In Sec. IIIII the experimental setup and the mea- 
surement procedures are described. Section IIVI presents 
various collective coupling rate measurements and com- 
pares them to the theoretical expectations. Section IVl 
shows measurements of the coherence time of collective 
coherences between Zeeman sublevels. A conclusion is 
given in Sec. IVII 



II. CQED INTERACTION: THEORETICAL 
BASIS 

A. Hamiltonian and evolution equations 




FIG. 1. Scheme considered for the description of the inter- 
action between an ion Coulomb crystals and the cavity field. 
The cavity is formed by the two mirrors Mi and M2 . eti n is the 
input light field, a the intracavity field, and a r is the reflected 
field, k denotes the cavity field decay rate. The spontaneous 
dipole decay rate of the ions is denoted by 7. 

We consider the interaction of -/V tot two-level ions in a 
Coulomb crystal with a single mode of the electromag- 



netic field of an optical cavity (denoted by nm), as de- 
picted in Fig. [T] The single-ended linear cavity is formed 
by two mirrors Mi (partial transmitter, PT) and M2 
(high reflector, HR) with intensity transmission coeffi- 
cients Ti and T2 {T\ ^> I2). The absorption loss co- 
efficient per round-trip is C and the empty cavity field 
round-trip time is r = 2//c, where I is the cavity length 
and c the speed of light. The intracavity, input and 
reflected fields are denoted by a, a- ln , and a r , respec- 
tively. The interaction of an ensemble of N identical 
two-level ions with a single mode of the cavity field can 
be described by a Jaynes-Cummings Hamiltonian of the 
form (HH 



H — H at + Hi + H a 



(1) 



where, in the frame rotating at the laser frequency uji, 
the atom and light Hamiltonians are given by H at = 
/lA^^j* fr^ and Hi = h/S. c a)a. The atomic and cavity 
detunings are denoted by A = cj at — wi and A c = uj c — u\, 
where aj at and w c are the atomic and cavity resonance 

~ (e) 

frequencies, respectively. 7rJ- is the excited state popula- 
tion operator of the j-th ion and a, are the intracavity 
field annihilation and creation operators. In the rotating 
wave approximation the interaction Hamiltonian reads 



H, 



al 



N t 



(2) 



where o\J and &j are the atomic rising and lowering oper- 
ators, defined in the frame rotating at the laser frequency. 
The single- ion coupling rate g is defined as g = /x ge -Eo / K 
where /i g0 is the dipole element of the transition consid- 
ered and Eq the maximum electric field amplitude. The 
field distribution -Eo^nm^j) is assumed to be that of a 
single-cavity Hermite-Gauss mode [7(j. In the follow- 
ing, we will restrict ourselves to the fundamental TEMoo 
mode of the cavity and refer to Ref. [56| for the coupling 
of ion Coulomb crystals to higher-order cavity transverse 
modes. 

The coupled atom-cavity system is subject to decoher- 
ence, mainly through the spontaneous decay of the ions 
from the excited state and through the decay of the cav- 
ity field due to the finite reflectivity of the cavity mirrors 
and due to intracavity losses. These dissipativc processes 
are characterized by the atomic dipole decay rate, 7, and 
by the total cavity field decay rate, k, respectively. The 
cavity field decay rate is given by k = n\ + «2 + «£, and 
includes the decay rates through the PT and HR mirrors 
(«i = Ti/2r and K2 = T2/2T) and the decay rate due to 
absorption losses (kc = £/2t). 

We derive standard semiclassical equations of motion 
for the mean values of the observables via (a) = ^([H, a]) 
and phenomenologically adding the relevant dissipativc 
processes @, i, H| MuM ■ 

In the low saturation regime, most of the atoms re- 
main in the ground state, (th ) <ti 1 , and the dynamical 
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FIG. 2. (Color online) (a) Calculated probe reflectivity spec- 
trum as a function of cavity detuning A c for an empty cavity 
(red line) and for N = 500 ions interacting with the cav- 
ity field (blue line) [see Eq. (jTHJ)] for a probe detuning of 
A = 7. The coupling of the ions and the cavity mode lead 
to a broadening and a shift of the resonance dip. The values 
of the parameters are set to typical values used in the exper- 
iments presented in Sees. IIIII and IIVI n = 2n x 2.1 MHz, 
«i = 2tt x 1.5 MHz, 7 = 2tt x 11.2 MHz, g = 2ir x 0.53 MHz. 
(b) Expected effective cavity decay rate «' and shift of the 
cavity resonance A^ — A c as a function of probe detuning A 
calculated for the same set of parameters. 



equations for the mean values of the observables read 
&j = -(7 + iA)<jj + ig^oo(rj)a, (3) 

JVtot 

a = -(« + iA c )a + i ^ g# o{rj)(Tj + v / 2Ki/ra in (4) 
where o = (6) is the mean value of observable 6. 



B. Steady-state reflectivity spectrum and effective 
number of ions 

In steady-state, the mean value of the intracavity field 
amplitude is given by 



iA> 



(5) 



where an effective cavity decay rate and an effective cav- 
ity detuning arc introduced: 



k' = k + g 2 N 



A' C = A C - g 2 N 



7 2 + A 2 ' 
A 



7 2 + A 2 



(6) 
(7) 



In these expressions, N is the effective number of ions 
interacting with the intracavity field, which is calculated 
by summing over all ions and weighting the contribution 
of each ion by the field modefunction under consideration 
evaluated at the ion's position: 



N 



iVtot 

E 

3=1 



Here, 



wo 
w(z) 

2 [k. 



exp 



*oote). 



2r 2 



(8) 



x sin 



w(z) 2 _ 

arctan(2/z ) + kr 2 /2R(z)] , (9) 



is the modefunction of the cavity fundamental TEMoo 
Gaussian mode with waist wq at the center of 
the mode and w(z) — Woyl + z 2 / ' z 2 , R(z) = z + z 2 /z, 
zq = 7tui 2 /A, and k = 2ir/\. 

Large ion Coulomb crystals in a linear radiofrequency 
trap are to an excellent approximation spheroids with 
half-length L and radius R (see Fig. [5]) , where the den- 
sity of ions, p, is constant throughout the crystal (53l.l54j]. 
It is then convenient to adopt a continuous medium de- 
scription, in which Eq. (j8]) becomes an integral over the 
crystal volume V: 



N = p f dr* 2 (r) 
Jv 



(10) 



In our experiment, the crystal radius and half-length, 
R and L, are typically much smaller than zq and the 
axial mode function can be approximated by sin 2 (fc2;). 
Moreover, for randomly distributed ions along the cavity 
axis z, one can average over the cavity standing- wave 
longitudinal structure, which gives an effective number 
of ions equal to 



N= - [ drexp[-2r 2 /w(z)' 2 
2 Jv 



(11) 



This expression can be evaluated knowing the crystal di- 
mensions, its density, and the cavity mode geometry. For 
typical crystals with large radial extension as compared 
to the cavity waist R 3> wo and length smaller than the 
Raylcigh range L <C zq, this expression reduces to 



N 



—L, 



(12) 



which is simply the product of the ion density by the 
volume of the cavity mode in the crystal. 
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Using the input-output relation a r = \J1k,\t<i— a- m , one 
finds that the steady-state probe reflectivity spectrum of 
the cavity is also Lorcntzian-shapcd in presence of the 
ions, the bare cavity decay rate, and detuning k and A c 
being replaced by their effective counterparts k' and A' c 
of Eqs. USE): 



a r 


i 







2«i - k' - iA' 



k' + iA' 



(13) 



The broadening and shift of the cavity resonance then 
represent the change in absorption and dispersion expe- 
rienced by the cavity field interacting with TV ions. In 
Fig. [2] (a) the expected cavity reflectivity spectrum is 
shown for both an empty cavity and a cavity containing 
a crystal with an effective number of ions N = 500 and 
for parameters corresponding to those used in the exper- 
iments presented in Sees. IIIII and IIVI In Fig. [2] (b) the 
effective cavity decay rate, k', and the shift of the cav- 

I 



ity resonance induced by the interaction with the ions, 
A' c — A c , are shown as a function of the probe detuning, 
A, for the same parameters. 



C. Effect of the motion of the ions 

The interaction Hamiltonian in Eq. ([2]) is only valid for 
atoms at rest. If an ion is moving along the axis of the 
cavity, the standing- wave structure of the cavity field and 
the Doppler shifts due to the finite velocity of the ion have 
to be taken into account. For an ion moving along the 
standing wave field with a velocity Vj , it is convenient to 
define atomic dipole operators, <jj± = iaj exp (±ikzj), 
arising from the interaction with the two counterpropa- 
gating components of the standing-wave cavity field. In 
the low saturation limit and taking into account the op- 
posite Doppler-shifts, the evolution equations ©,(|3} be- 
come 



< J j± 



[7 + i(A ± kvj)} <r j± + i(g/2)y nm ( rj )a 



N t , 



a = -(« + iA c )a + i(g/2) ^ *nm(rj) (<r j+ + + y / 2K 1 /Ta i: 



(14) 
(15) 



When the typical timescale of the motion is slow as com- 
pared to the timescales for the coupled dynamics of the 
atomic dipole and cavity field, the steady-state mean 
value of the intracavity field can be found by averaging 
the contributions of the individual dipole mean values 
given by Eq. (|14[) over the distribution of the mean ve- 
locities, f(v). For a distribution f(v) with an average 
velocity ud a conservative estimate for this to be valid is 
that the mean Doppler-shift is smaller than both effective 
rates of the coupled system on resonance (A c = A = 0), 
kvn <C min[ft + g 2 N/~f,"f + g 2 N/K}. Under these condi- 
tions, the expression for the intracavity field mean value 
is then of the same form as in the zero-velocity case 
(Eq. ([5])). The effective cavity field decay rate and de- 
tuning of Eqs. ([H]) and© are modified according to 



k' = k + g 2 N / dvf(v)~/£(v) 



A' C = A C - g 2 N J Avf{v){A - kv)£(v). (17) 



where 



1 2 + A 2 + {kvf 



( 7 2 + A 2 ) 2 + 2( 7 2 - A 2 ){kv) 2 + (kv) 4 



(18) 



In the case of a thermal Maxwell-Boltzmann 

distribution with temperature T, one has 

= v/5^T CX P whcre fc B is the Boltz- 

mann constant and m the mass of the ion. At low 



temperatures, i.e., when the width of the thermal 
distribution is small as compared to the atomic natural 
linewidth, the effective cavity field decay rate and detun- 
ing given by Eqs. (fll)j) and (|T7|) are well-approximated 
by 



/ 1 2 AT 7 

y 7 ' 2 + A 2 ' 
A^ = A c - g 2 N- 



(19) 
(20) 



7 ' 2 + A 2 ' 

These equations are of the same form as Eqs. (|6|) and (|7|), 
replacing the natural dipole decay rate by an effective 
dipole decay rate , 



7'~7(l + fcv D /V2), 



(21) 



(16) 

where vr, = \JksT jm is the mean Doppler velocity. 



III. EXPERIMENTAL SETUP 

A. Cavity trap 

The ion trap used is a segmented linear quadrupolc 
radiofrequency trap that consists of four cylindrical elec- 
trode rods (for details see [zlj]). The electrode radius is 
2.60 mm and the distance from the trap center to the 
electrodes is r = 2.35 mm. Each electrode rod is di- 
vided into three parts, where the length of the center 
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FIG. 3. (Color online) Schematic experimental setup. The ab- 
breviations are: polarizing beam splitter (PBS), single mode 
fiber (SMF), acousto-optical modulator (AOM), dichroic mir- 
ror (DM), Pound-Drever-Hall lock (PDHL), avalanche photo- 
diode (APD), second harmonic generation (SHG). The pho- 
toionization laser is not shown. 



electrode is zq — 5.0 mm, and the length of the end elec- 
trodes is ze = 5.9 mm. Radial confinement is achieved 
by a radiofrequency field (RF) applied to the entire rods 
at a frequency of 2ir x 4 MHz and a it phase difference 
between neighboring rods. The axial trapping potential 
is created by static voltages (DC) applied to the outer 
parts of the rods. 

An optical cavity is incorporated into the trap with its 
axis parallel to the symmetry axis of the ion trap (see 
Fig. [3]). The cavity mirrors have a diameter of 1.2 mm 
and a radius of curvature of 10 mm. The rear face of 
both mirrors are anti-reflection coated at a wavelength of 
866 nm corresponding to the 3d 2 D 3 / 2 <-> 4p 2 P 1 / 2 transi- 
tion in 40 Ca + , while the front facade of one mirror is par- 
tially transmitting (PT) and for the other highly reflect- 
ing (HR) at this wavelength. Their intensity transmission 
coefficients are 1500 and 5 ppm, respectively. The intra- 
cavity losses due to contamination of the mirrors during 
the initial bake out amount to ~ 650 ppm. The PT mir- 
ror is mounted on a plate that can be translated using 
piezoelectric actuators to allow for scanning or actively 
stabilizing the cavity length. The cavity has a close to 
confocal geometry with a length of 11.8 mm, correspond- 
ing to a free spectral range of 12.7 GHz and a waist of the 
fundamental TEMoo mode of wq = 37/im. With a mea- 
sured cavity field decay rate of k = 2tt x (2.1 ± 0.1) MHz, 
the finesse is found to be T = 3000 ± 200 at a wavelength 



of 866 nm [74|. 

40 Ca + ions are loaded into the trap by in situ photoion- 
ization of atoms from a beam of atomic calcium in a two- 
photon resonant photoionization process [74|-[76j. The 
ions are cooled to a crystalline state through Dopplcr- 



laser cooling using a combination of two counterpropa- 
gating laser beams, resonant with the 4s 2 S!/ 2 <-> 4p 2 Pi/ 2 
transition at 397 nm along the trap axis, and a repump- 
ing laser applied along the x axis and resonant with 
the 3d 2 D 3/ / 2 4p 2 P!/ 2 transition at 866 nm to pre- 
vent shelving to the metastable D 3 / 2 state. Three sets 
of Helmholtz coils are used to compensate for residual 
magnetic fields and to produce bias magnetic fields. For 
the measurements of the collective coupling rate between 
the ion Coulomb crystals and the cavity light field, the 
transverse magnetic fields along x and y are nulled and 
a magnetic field of B z ~ 2.5 G along the z-axis is used. 



B. Detection 

A grating stabilized diode laser at 866 nm provides the 
light for probing the coupling of the ion Coulomb crystals 
with the standing wave field inside the optical cavity. It 
is injected into the cavity through the PT mirror. Ad- 
ditionally, a second grating stabilized diode laser with a 
wavelength of 894 nm serves as an off-resonant reference 
laser and is simultaneously coupled to the cavity through 
the PT mirror and used to monitor the cavity resonance. 
Both lasers are frequency stabilized to the same temper- 
ature stabilized reference cavity and have linewidths of 
- 100 kHz. 

The reflectivity of the 866-nm cavity field is measured 
using an avalanche photo diode (APD). The light sent 
to the APD is spectrally filtered by a diffraction grat- 
ing (1800 lines/mm) and coupled to a single mode fiber. 
Taking into account the efficiency of the APD at 866 nm, 
the fiber incoupling and the optical losses, the total de- 
tection efficiency amounts to ~ 16%. A similar detection 
system is used to measure the transmission of the 894- nm 
reference laser. 

Depending on the experiment, the reference laser 
serves two different purposes. In a first configuration, the 
length of the cavity is scanned at a rate of 30 Hz over the 
atomic resonance. In this configuration, the frequency 
of the reference laser is tuned such that it is resonant at 
the same time as the probe laser in the cavity scan. This 
allows for monitoring slow drifts and acoustic vibrations. 
The signal of the weak probe laser is then averaged over 
typically 100 scans in which the stronger reference laser 
is used to keep track of the current position of the cavity 
resonance. 

In a second configuration, the cavity resonance is 
locked on the atomic resonance by stabilizing the length 
of the cavity to the frequency of the reference laser in a 
Pound-Drever-Hall locking scheme [77j • During the mea- 
surement, imperfections in the stabilization are compen- 
sated for by monitoring the transmission of the 894-nm 
reference laser. The data is then postselected by only 
keeping data points for which the transmitted reference 
signal was above a certain threshold. 
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C. Experimental sequence 

In both configurations, the cavity reflection spectrum 
is measured at a rate of 50 kHz using a 20 fj,s sequence 
of Doppler cooling, optical pumping and probing, as in- 
dicated in Fig. 2J First, the ions are Doppler- laser cooled 
for 5 fis by driving the 4s 2 S!/ 2 <-> 4p 2 P 1 / 2 transition us- 
ing laser cooling beams at 397 nm (LC), and at the same 
time repumping on the 3d 2 D 3 / 2 <-> 4p 2 P]y 2 transition 
with a laser at 866 nm (RP). Next, the ions are opti- 
cally pumped to the mj = +3/2 magnetic substate of 
the 3d 2 D 3 / 2 level by applying the optical pumping laser 
(OP) in combination with the laser cooling beams (LC) 
for a period of 12 ps. The optical pumping laser is res- 
onant with the 3d 2 D 3 / 2 <-> ^P 2 Pi/2 transition and has a 
polarization consisting only of <7 + - and 7r-polarized com- 
ponents. It is sent to the trap under an angle of 45° with 
respect to the quantization axis. By probing the popula- 
tions of the different Zeeman sublevels, the efficiency of 
the optical pumping was measured to be rj = 97lg% [7a [. 
Finally, the cavity reflection signal is probed by inject- 
ing a 1.4 ps cr~-polarized probe pulse, resonant with the 
3d 2 D 3 / 2 <-> 4p 2 P!/ 2 transition, into the TEM o mode 
of the optical cavity. Its intensity is set such that the 
mean intracavity photon number is less than one at any 
time. With a delay of 0.1 ps relative to the probe laser, 
the APD is turned on. The delay ensures that the field 
has built up inside the cavity and that the system has 
reached a quasi-steady-state. The length of the probing 
period was chosen in order to minimize the total sequence 
length as well as to avoid depopulation due to saturation 
of the transition [7811. 



D. Effective number of ions 

As mentioned above, the effective number of ions in- 
teracting with the cavity field depends on the ion crystal 
density and the overlap between the crystal and the cav- 
ity modevolumc, where the density of the ion Coulomb 
crystals depends on the amplitude of the RF voltage [55[ : 



lated using the formula 

N = r\t J dxdy cxp {-2[(x - x ) 2 + (y - y ) 2 ]/w 2 } , 

(23) 

where the parameter r/ accounts for a finite efficiency of 
the optical pumping preparation and xq and yo denote 
the radial offsets. These offsets can in principle be can- 
celed to within a micron [80j . but in the experiments re- 
ported here, they were measured to be xq = 3.9 /im, yo = 
15.7 /im [56j. The uncertainty in the effective number of 
ions comes from the uncertainty Sp in the density de- 
termination, due to the RF voltage calibration, the un- 
certainty in the crystal volume 6V, due to the imaging 
resolution Sx and the uncertainty of the optical pump- 
ing efficiency Srj. The relative uncertainty in the effective 
number of ions, TV = rjpV, can then be expressed as (78j 




SV 
~~V 



(24) 



where 6V/V = 6x\/16L 2 + R 2 /2RL. For the typically 
few-mm-long prolate crystals used in these experiments 
and an imaging resolution 5x ~ pm, this results in a 
relative uncertainty of 5-7% in the effective number of 



IV. COLLECTIVE COUPLING 
MEASUREMENTS 

To achieve collective strong coupling on the chosen 
3d 2 D 3/2 , mj = +3/2 O 4p 2 P 1/2 , mj = +1/2 transi- 
tion the collective coupling rate gy/N has to be larger 
than the cavity field decay rate k — 2ir x 2.1 MHz and 
the optical dipolc decay rate 7 = 2n x 11.2 MHz. With 
the known dipole element of the transition and the cavity 
geometry, the single-ion coupling rate at an antinode at 
the center of the cavity fundamental mode is expected to 
be g = 2?r x (0.53 ± 0.01) MHz. One thus expects to be 
able to operate in the collective strong coupling regime 
as soon as TV > 500. 



£ o^rf 
1U 'o"rf 



(22) 



A. Atomic absorption and dispersion 



Here, M denotes the ion mass. The precise calibration of 
the RF voltage on the trap electrodes can be performed, 
e.g., on the basis of a zero-temperature charged liquid 
model [HI, [79|. [80l or the measurement of the Wigner- 



Seitz radius [801 ]. For the trap used in these experiments, 
p = (6.01 + 0.08) x 10 3 U&p V~ 2 cm- 3 . The crystal mode 
volume is found by taking fluorescence images of the crys- 
tal during Doppler-laser cooling, as shown in Fig. [SI from 
which the crystal half-length L and radius R can be ex- 
tracted. Taking a possible offset between the cavity axis 
and the crystal revolution axis into account, the effective 
number of ions [see Eq. (|lll) ] is then numerically calcu- 



To investigate the coherent coupling of the ions with 
the cavity field in the collective strong coupling regime, 
we first perform measurements of the atomic absorp- 
tion and dispersion of a given crystal with TV ~ 500 by 
scanning the cavity length around atomic resonance and 
recording the probe reflectivity spectrum. The crystal 
used in these experiments is similar to the one shown in 
Fig. El With a density of p = (5.4 ± 0.1) x 10 8 cm- 3 
, a half-length L = (511 + 1) pm and radius R = 
(75 ± 1) pm the total number of ions in the crystal is 
Ntot = 6500 ± 200, and the effective number of ions in- 
teracting with the cavity mode is TV = 520 i 3 2. 
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FIG. 4. (color online) (a) Experimental sequence used to measure the collective coupling rate, (b) Energy levels of 40 Ca + 
including the relevant transitions addressed in the three parts of the experimental sequence. The acronyms are: laser cooling 
beam (LC), repumping beam (RP), optical pumping beam (OP), probe beam (PB), avalanche photodiode (APD). 




FIG. 5. (Color online) (a) Typical projection image of a 
crystal used in collective strong coupling measurements. All 
ions are exposed to cooling and repumping light. The solid 
blue line indicates the outline of the ellipsoid, (b) Same 
crystal, but only the ions in the cavity mode are exposed 
to repumping light, now injected into the cavity. The ions 
outside the cavity mode volume are not visible, as they are 
shelved into the metastable 3d 2 D 3 / 2 level. The crystal con- 
tains iVtot = 8780 ± 180 ions, of which N = 4891^ effectively 
interact with the cavity field. 



The broadening and the shift of the cavity resonance 
are then measured as a function of the detuning of the 
probe laser, A. This is accomplished by scanning the cav- 
ity length over a range corresponding to ~ 1.3 GHz at a 
repetition rate of 30 Hz, for a fixed value of A. The width 
of the reflection dip for a given detuning A is found by 



averaging over 100 cavity scans, where the reference laser 
is overlapped with the probe laser on the cavity scan and 
used to compensate for any drift of the cavity. In Fig. [51 
cavity reflection scans are plotted for various detunings. 
Each data point corresponds to the average of 100 20— /is- 
measurement sequences as showed in Fig. HI As expected 
from Eq. (|19p . the broadening of the intracavity field ab- 
sorption reflects the two-level atomic medium absorption. 
Each set of data is, according to Eq. (fT3|) . fitted to a 
Lorentzian from which the cavity half width half max 
imum (HWHM) k 
modified cavity HWHM, k 
of the probe laser, A. Each point is the average of five 
measurements; the solid line is a fit according to Eq. (|19l) . 
From the fit we deduce a collective coupling rate of <7n = 
27r x (12.2 ± 0.2) MHz, in good agreement with the the- 
oretical expectation of <7n, theory 

= 2tt x (12.1±g;f) MHz, 
calculated for N = 520^32 ions interacting with the cav- 
ity mode [28|]. Furthermore, the effective dipole decay 
rate 7' is left as a fit parameter to account for nonzero 
temperature effects, as discussed in Sec. Ill CI The fit 
yields 7' = 27r x (11.9 ± 0.4) MHz, which would corre- 
spond to a temperature of T = 24+jjj mK, and a natural 
half-width of the cavity of K = 2n x (2.2 ± 0.1) MHz, in 
good agreement with the value deduced from an indepen- 
dent measurement of the free spectral range (FSR) and 
the finesse of the cavity, k = 2ttx (2.1 ± 0.1) MHz 

For the measurement of the effective cavity detuning, 
A' c , the position of the 894- nm resonance laser in the 
cavity scan is fixed to the bare cavity resonance. The 
frequency shift is then measured by comparing the posi- 



' is deduced. Figure [Tfa) shows the 
as a function of detuning 
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FIG. 6. (Color online) Typical probe reflectivity for various 
values of the atomic detuning A. The probe detunings were 
(a) A « 2tt x 54.3 MHz, (b) A « 2tt x 24.3 MHz, (c.) A « 
2tt x 8.3 MHz, and (d.) A^2ix 0.3 MHz. Solid lines are 
Lorentzian fits to the data; the effective cavity field decay rate 
k is deduced from the fit. 
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FIG. 7. (Color online) (a) Measured cavity field effective 
decay rate ft' versus probe detuning A for a crystal with 
N = 520t 2 3 t ions interacting with the cavity field. The blue 
solid line is a fit to the data, (b) Measured shift of the cavity 
resonance frequency A^ — A c versus atomic detuning A for 
the same crystal. The blue line is a fit to the data. 



tion of the probe and the reference signal resonances in 
the cavity scan. The effective cavity detuning as a func- 
tion of probe detuning is shown in Fig.JTJb) One observes 
the typical dispersive frequency-shift of two-level atoms 
probed in the low saturation regime. The data is fitted 
to the theoretical model according to Eq. (|20|) . to find a 
collective coupling rate <7n = 27r x (12.0 ± 0.3) MHz and 
an effective dipole decay rate 7' = 27r x (12.7 ±0.8) MHz. 
Both values are consistent with the previous measure- 
ment and the theoretical expectations. As in the pre- 
vious measurement, the 894-nm reference laser is used 
to compensate systematic drifts and acoustic vibrations. 
However, since this compensation method relies on the 
temporal correlations of the drifts in both signals, and 
thereby on their relative positions in the cavity scan, the 
compensation becomes less effective at large detunings. 
This is reflected in the bigger spread and the larger er- 
ror bars at larger detunings, which renders this method 
slightly less precise than the absorption measurement to 
evaluate the collective coupling rate. 
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Probe and cavity detuning A = Ac (2ir MHz) 



FIG. 8. (color online). Probe reflectivity signal as a function 
of A = A c for the empty cavity (blue triangles) and with a 
crystal with N = 520^32 effectively interacting ions present 
in the cavity mode volume (red circles). The solid lines are 
fits to the theory. 



B. Vacuum Rabi splitting 

A third complementary method to measure the 
collective coupling rate is based on locking the cavity 
on atomic resonance, lu c = w a t. The response of the 
coupled atom-cavity system is then probed as a function 
of probe detuning A, which is then equal to the cavity 
detuning A c . The result of this measurement is shown 
on Fig. [5] The blue triangles are obtained with an 
empty cavity, while the red circles were taken with 
the same ion Coulomb crystal as used in the previous 
experiments. Each data point is deduced from 2 x 10 4 
experimental sequences (see Fig. [4j. The results are 
fitted using the theoretical expectations of Eq. (fT3|) 
and Eqs. ([19) and ([20]) (solid lines in Fig. H) and yield 
<7n = 27T x (12.2 ± 0.2) MHz, a value that is in good 
agreement with the previous measurements. To facilitate 
the convergence of the more complex fitting function, 
the value of 7' in Eqs. ([T9| and ([20]) was set to the one 
found in the previous absorption measurement. 

From these three independent measurements of the 
collective coupling rate and using the effective num- 
ber of ions one deduces a single ion cou- 
pling rate of g cxp = 2n x (0.53 ± 0.02) MHz, which 
is in excellent agreement with the expected value of 
5thcory = 2tt x (0.53 ± 0.01) MHz. 



C. Scaling with the number of interacting ions 

To check further the agreement between the theoret- 
ical predictions and the experimental data, we investi- 
gated the dependence of the collective coupling rate on 




Effective number of ions N 



FIG. 9. (Color online) Cooperativity as a function of the ef- 
fective number of ions. The solid line is a linear fit to the data 
and yields a scaling parameter of = (5-IIq' 4 .) x 10 -3 . The 
shaded area indicates the confidence region for the upper and 
lower limit of C determined by the (systematic) uncertainties 
on the effective number of ions, N. The dashed line indicates 
the strong collective coupling limit <?n > (^,7). 



the effective number of ions. An attractive feature of ion 
Coulomb crystals is that the number of ions effectively 
interacting with a single mode of the optical cavity can 
be precisely controlled by the trapping potentials. While 
the density p only depends on the amplitude of the RF 
voltage [sec Eq. (|22|)]. the aspect ratio of the crystal de- 
pends on the relative trap depths of the axial and radial 
confinement potentials, which can be independently con- 
trolled by the DC voltages on the endcap electrodes. This 
allows for controlling the number of effectively interact- 
ing ions down to the few ion-level. 

By analogy with the case of a single two-level sys- 
tem interacting with a single field mode of an optical 
cavity, the cooperativity parameter C is defined here as 
(half) the ratio of the square of the effective coupling 
rate <7n to the cavity field decay rate k times the effec- 
tive dipole decay rate 7' (taking into account the effect 
of the motion of the ions): C = g^/2n"f'. As can be seen 
from Eq. (|19[) . this parameter can be experimentally ob- 
tained by measuring for a probe field tuned to atomic 
resonance (A = 0) the effective cavity field decay rate 

k'(A = 0) = k+ ^- = k(1 + 2C). In Fig. H the de- 
pendence of the cooperativity parameter, C, is plotted 
as a function of the effective number of ions interacting 
with the TEMoo mode, where the effective number of ions 
was changed by measuring for different aspect ratios and 
densities of several crystals. 

The effective number of ions in each crystals was de- 
duced by applying the method described in Sec. IIIIDI 
The data points were obtained using er~-circularly po- 
larized probe light, hence probing the population in the 
toj = +3/2 and toj = +1/2 substates, and shows the 
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FIG. 10. (color online). Vacuum Rabi splitting spectra (A = 
A c ) obtained for increasing effective number of ions [0 (blue 
squares), 243 (lilac circles), 601 (dark red triangles), 914 (red 
diamonds)], the lines are presented to guide the eye. 

expected linear dependence on the effective number of 
ions. From a linear fit (solid line) we deduce a scaling 
parameter = (5.1^02) x 10~ 3 . The limit where collec- 
tive strong coupling is achieved (<?n > k, 7) is indicated 
by the black dashed line and is reached for ss 500 inter- 
acting ions. 

The largest coupling observed in these experiments was 
measured for a crystal with a length of ~ 3 mm and 
a density of ~ 6 x 10 8 cm" 3 and amounted to C = 
7.9 ± 0.3, corresponding to an effective number of ions of 
N = 1 523!gg . This value exceeds previously measured 
cooperativities with ions in optical cavities by roughly 
one order of magnitude [14l - tl6l |. 

Similarly, vacuum Rabi splitting spectra, such as the 
one presented in Fig. [8j were measured for several crys- 
tals and aspect ratios. The result of such measurements 
is shown in Fig. 1101 showing clearly the increase in the 
separation between the coupled crystal+cavity normal 
modes as the number of ions is increased. The collec- 
tive coupling rate g-^, derived from fits to the theoretical 
expression Eq. (|13[) . is plotted for different effective num- 
ber of ions in Fig. 111! Taking the finite optical pumping 
efficiency into account and fitting the curve with the ex- 
pected square-root dependency, we deduce a single ion 
coupling rate of g = 2tt x (0.53 ± 0.01) MHz, in good 
agreement with the previous measurements and the the- 
oretical expectation. 




200 400 600 800 1000 
Effective number of ions N 



FIG. 11. (color online). Collective coupling rate versus 
effective number of ions N deduced from reflectivity spectra, 
such as shown in Fig. 1101 obtained with crystals of different 
shape and density. The blue line is a fit to the data and 
gives a single ion coupling rate g = 2tt x (0.53 ± 0.01) MHz. 
The shaded area indicates the lower and upper bound of the 
collective coupling rate g^ within the uncertainties of N. The 
horizontal errorbars are calculated according to Eq. (|24l) 

level. These coherences were established by the Larmor 
precession of the magnetic spin induced by an additional 
_B-ficld transverse to the quantization axis. In presence of 
this orthogonal i?-field, the population of the several sub- 
states undergo coherent oscillations, which are measured 
at different times in their free evolution by directly prob- 
ing the coherent coupling between the cavity field and 
the ions. In order to be able to resolve the coherent pop- 
ulation oscillations in time using the previous technique 
(probing time ~ 1 /is) the amplitude of the longitudi- 
nal -B-field was lowered to obtain oscillation periods in 
the ~ 10 /is range, and the optical pumping preparation 
was modified as to minimize the effect of the transverse 
.B-field. The reduced .B-field along the quantization axis 
could in principle make the sample more sensitive to B- 
ficld fluctuations. Since these fluctuations might be one 
of the factors eventually limiting the achievable coher- 
ence time, we expect the coherence time measured by 
this method to be a lower bound as compared to the 
previous configuration with a larger longitudinal B-ficld. 



A. Experimental sequence and theoretical 
expectations 



V. COHERENCE TIME OF COLLECTIVE 
ZEEMAN SUBSTATE COHERENCES 

To evaluate the prospect for realizing coherent manip- 
ulations, we measured the decay time of the collective 
coherences between the Zeeman substates of the 3d 2 D 3 / 2 



The coherence time measurements required the exper- 
imental configuration and the measurement sequence to 
be slightly modified as compared to the collective cou- 
pling rate measurements described in Sec. IIII CI The 
Larmor precession is induced by an additional .B-field 
component along the transverse x direction, while the 
longitudinal magnetic field component B z was lowered 
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FIG. 12. (Color online) (a) Experimental sequence used to measure the coherence time of collective Zeeman substate coherences 
in the 3d 2 D 3 / 2 level, (b) Energy levels of 40 Ca + , including the relevant transitions and their polarization for the four phases 
of the experimental cycle. In the third phase, all lasers are turned off for a variable delay r, and the system evolves freely in 
presence of a transverse magnetic field component _B X . 



to optimize the contrast of the coherent population oscil- 
lations. The optical pumping light propagates along the 
x axis and is 7r-polarized, hence transferring most of the 
atoms symmetrically into the two outermost magnetic 
sub-states of the 3d 2 D 3 / 2 level, mj = ±3/2. 

The experimental sequence used to measure the coher- 
ence time is shown in Fig. [12] The ions are Doppler-laser 
cooled during the first 5 /is, followed by a 12 /is optical 
pumping period. After the optical pumping, all lasers are 
turned off for a time t, allowing for the free evolution of 
the system. Finally, a weak cr _ -circularly polarized probe 
pulse is injected into the cavity, addressing the ions in the 
rnj = +1/2 and mj = +3/2 sub-states. The steady-state 
cavity reflection is measured by collecting the reflected 
photons with the APD for 0.5 fis. The additional delay 
time between optical pumping preparation and probing 
obviously lowers the repetition rate of the sequence sig- 
nificantly, especially for long delay times, and the number 
of data points for each sweep of the cavity will decrease. 
To compensate for this, the data points at longer delays 
had to be averaged over more cavity scans, which sub- 
stantially increased the acquisition time and eventually 
limited these measurements to delays of around ~ 120 /is. 

Based on a simple four-level model the free Larmor 
precession-induced changes in the populations of the Zee- 
man substates, |mj = ±1/2, ± 3 /2), of the 3d 2 D 3 / 2 level 
can be calculated. For a homogeneous £>-ficld with com- 
ponents -B x and -B z , the Hamiltonian of the four- level 
system can be expressed in terms of collective popula- 



tion operator, 

6^ = ^^^/', (25) 
and collective spin operators 

N tot 
J'=l 

Here, \mj)^' and \m! } )^ are the state kets of the jth ion 
with magnetic quantum number m j and m' } , respectively. 
The sum extends over the total number of ions. In this 
notation, the Hamiltonian of the free evolution of a spin 
J = 3 /2 system can be written as 

H B = hw z ^2 m jn m ., + x 

m .l( TO J - l)^mj,rai + l+ (27) 

V T ~ mj ( mj + 1 ) < Wmj- 1 ) 

where the sums extend over the four Zeeman substates. 
Here, 5 mj>m '_|_i is the Kronecker delta, and the Larmor 
frequencies uj z and w x corresponding to the z and x com- 
ponent of the magnetic field arc given by the product of 
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the magnetic field amplitude by the gyromagnetic ratio 
7gm: 

w z = 7gm-B z , w x = 7gm-B x . (28) 

For a cr~-circularly polarized probe, the measured collec- 
tive coupling to the cavity light will depend on the col- 
lective populations in the to,j = +1/2 and mj = +3/2 sub- 
states. For a nonvanishing population in the mj = + l /2 
state, the measured effective cavity decay rate, which was 
defined for a two-level system in Eq. dH]), contains both 
contributions and is hence modified to 

«'(t) =W + g? /a iVy a (T) 7 A2 +ff3/2%2(T) 2 7 A2 . 

' V 2 ' 3 / 2 

(29) 

where g> OTj , iV mj , and A mj = w mj — wi denote the single- 
ion coupling rate, the effective number of ions and the 
atomic detunings of the relevant Zeeman substates mj = 
+Y2, + 3 /2, respectively, and uj mj is the frequency of the 
3d 2 D 3 / 2 , raj -B- 4p 2 P!/ 2 ,raj — 1 transition. 

Due to the induced Larmor precession, the effective 
number of ions in the individual Zeeman substates will be 
time-dependent. For a system initially prepared in a su- 
perposition state i/}q, the population in a particular Zee- 
man substate at a certain time r can be calculated from 
the projection of the time evolved state, tjj( T ) = U(t)i/jo, 
onto this state. Here, U(t) = exp^t/nHsr) denotes 
the time evolution operator. Straightforward but lengthy 
calculations show that the populations in the + 1 /2 and 
+ 3 /2 Zeeman substates after a time r are of the form 
Acos{u>i,t)+B cos(2u;lt)-)-C, where A, B, and C are con- 
stants depending on the efficiency of the optical pumping 
(i.e., the initial populations and coherences in the differ- 
ent Zeeman sublevels) and the magnetic field amplitudes 
Bz and B x (via w x and uj z ). One thus obtains Aq/ 2 ( T ) 
and N 3 / 2 {t) using Eq. (8). It follows from Eq. ([2"Tf and 
k'(t) = k(1 + 2C(t)) that the measured cooperativity at 
time t can be put under the form 

C(t) = (zcos(u;lt) + 6cos(2ojl''") + c, (30) 

where the Larmor frequency 

UJ-L = yjul+ul (31) 

was defined. The parameters a, &, c are constants de- 
pending on the efficiency of the optical pumping prepa- 
ration, and the magnetic field amplitudes B z and B x . 

B. Experimental results 

The amplitudes of the magnetic fields, B K and B z , 
at the position of the ion crystal were calibrated by 
measuring the dependence of the Larmor frequency 
«l with the intensity of the current used to drive the 
transverse magnetic field coils [sec Eqs. (|2"g)l and (pTTj) ]. 
The obtained coupling as a function of t is shown 
for different currents 7 X on Fig. 1131 The curves are 
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FIG. 13. (Color online) Calibration of the Larmor frequency 
for different currents of the B x coils. Shown is the cooper- 
ativity as a function of delay time r for different transverse 
B-fields: J x = 10 mA (blue squares), J x = 16 mA (lilac stars), 
7 X = 26 mA (dark red circles), and J x = 36 mA (red trian- 
gles). The solid lines are fits according to Eq. (|30[) . 




I x (mA) 

FIG. 14. (Color online) Larmor frequency as a function of 
current through the -B x coils. The solid line is a fit of the 
form cjl = V^i + d 2 Ii an d we deduce o; z = 2n x (0.150 ± 
0.002) MHz and w x = 2tt x (5.5 ± 0.1) |3f x 7 X . 

fitted according to Eq. (|30j). yielding the individual 
Larmor frequencies. These frequencies are shown as a 
function of the current through the J5 X coils in Fig. [TJ] 
Using the gyromagnetic ratio 7gm = ^' B ^ 3/2 (/^b is the 
Bohr magneton, g 3 / 2 the Lande factor of the 3d 2 D 3 / 2 
level), we deduce the magnetic fields along the two axis 
B z = (0.134 ± 0.002) G and B x = (4.91 ± 0.09) f x 7 X . 

To achieve a large contrast, the measurement was car- 
ried out with moderate -B-field values B x = B z = 0.15 G 
and the variation of the cooperativity was measured for 
120 fis. To compensate for slow drifts during the mea- 
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assuming a Gaussian decay exp(— t2 /V|) m Eq. (|30|) yields 
a coherence time of r g = 0.5^q'2 ms - Due to the limita- 
tion of our measurement to time delays of r < 120 fis, it 
is at present not possible to distinguish between the two 
decay mechanisms. 

For comparison, the cooperativity as a function of 
probe delay, C(t), was measured with only the bias field 
along the quantization axis present (B x = 0, B z = 
0.15G), as shown in Fig. Q1)] b. Here, the values are 
normalized to the mean cooperativity averaged over all 
points (C). Within the error bars, the deduced coopera- 
tivities agree with a constant value of (C) = 1.43 ± 0.02 
(solid line). 

VI. CONCLUSION 



FIG. 15. (Color online) (a) Normalized cooperativity 
paramter as a function of delay r. Due to the presence of 
a non-zero B-field component orthogonal to the quantization 
axis (B z = _B X = 0.15 G), coherent L armor precessions are ob- 
served. Long term drifts are compensated by normalizing to 
the mean of one oscillation period. The solid line corresponds 
to a fit, assuming an exponential decay and yields a coherence 
time of r c = 1.7l o 8 ms. (b) Cooperativity as a function of 
delay with a B-field present only along the quantization axis 
(£? x = By = 0, B z = 0.15 G). The data points are normalized 
to the mean cooperativity of (C) = 1.43 ± 0.02. 



surement, each data point was normalized to the mean 
cooperativity, C, averaged over one oscillation period. 
The normalized cooperativity is shown in Fig. I15f a). to- 
gether with a fit of the form of ([3TJ)) , where decoherence 
processes are taken into account by multiplying the oscil- 
lating terms with an exponential decay term exp(— T /-r ), 
which would be expected, e.g., for a homogeneous broad- 
ening of the energy levels. From this fit, we deduce a 
coherence time of r c = 1.7l° ° 8 ms - This value is compa- 
rable to previously measured coherence times for single 
ions in linear Paul trap in equivalent magnetic field sen- 
sitive states [8lj and might be further improved by an 
active control of stray magnetic fields or state configu- 
rations that are less magnetic field sensitive. For inho- 
mogeneous broadening, due to magnetic field gradient 
over the crystal, the decoherence process would be bet- 
ter described by a Gaussian decay (HJ. Fitting the data 



To conclude, we have presented a detailed theoreti- 
cal and experimental analysis of the experiments of [28| , 
which demonstrated the possibility of using large ion 
Coulomb crystals positioned in a moderately high-finesse 
optical cavity to enter the collective strong-coupling 
regime of CQED. The excellent agreement between the 
experimental results including those of Ref. [28[ and 
the theoretical predictions, makes ion Coulomb crystals 
promising candidates for the realization of quantum in- 
formation processing devices such as quantum memo- 
ries and repeaters |29|. IHoll. Using, for instance, cavity 
EIT-based protocols [43. l83l-l86j. the obtained coupling 
strengths and coherence times could open up for the re- 
alization of both high-efficiency and long life-time quan- 
Moreover, the nice properties of ion 



57] 



turn memories 
Coulomb crystals also allow for the manipulation of com- 
plex multimode photonic informatio n [ 5 71 by exploiting 
the crystal spatial [56| or motional [63| degrees of free- 
dom. Ion Coulomb crystals in optical cavities have also 
great potential for the investigation of cavity optomc- 

and the observation of novel 
68i with cold, solid-like objects. 



chanical phenomena J6^ 
phase transitions [U [65| 
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